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ABSTRACT 

We report the discovery of HD 156668 b, an extrasolar planet with a minimum mass of 
Mp sini = 4.15 Mq. This planet was discovered through Keplerian modeling of precise radial ve- 
locities from Keck-HIRES and is the second super-Earth to emerge from the NASA-UC Eta-Earth 
Survey. The best-fit orbit is consistent with circular and has a period of P = 4.6455 d. The Doppler 
semi-amplitude of this planet, _R' = 1.89ms~^, is among the lowest ever detected, on par with the 
detection of GJ581e using HARPS. A longer period (Pw2.3yr), low-amplitude signal of unknown 
origin was also detected in the radial velocities and was filtered out of the data while fitting the 
short-period planet. Additional data are required to determine if the long-period signal is due to 
a second planet, stellar activity, or another source. Photometric observations using the Automated 
Photometric Telescopes at Fairborn Observatory show that HD 156668 (an old, quiet K3 dwarf) is 
photometrically constant over the radial velocity period to O.lmmag, supporting the existence of the 
planet. No transits were detected down to a photometric limit of ~3mmag, ruling out transiting 
planets dominated by extremely bloated atmospheres, but not precluding a transiting solid/liquid 
planet with a modest atmosphere. 

Subject headings: planetary systems — stars: individual (HD 156668) — techniques: radial velocity 



1. INTRODUCTION 

The search for low-mass planets is driven by a de- 
sire to observationally study the full range of plane- 
tary systems in order to better understand their for- 
mation, dynamics, composition, and diversity. We 
also seek Earth-like worlds of terrestrial composition 
as a goal in itself and as targets for future transit 
and imaging observations. This search has taken sev- 
eral leaps forward recently because of instrumental im- 
provements. The precision of radial velocity (RV) mea- 
surements with Keck- HIRES by the California Planet 
Survey (CPS) gr oup (IHoward et al . 2009b ), HARPS 
(|Loviset al.ll2006l ). and the AAT ([O'Toole etan[200l 
has now reached lms~^ or better and has led to the 
discovery of several super-Earths around nearby, bright 
stars. Ground-based tran sit surveys such as MEarth 
(jCharbonneau et all |2009| 1 and HATNet (jBakos et al.l 
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l2009f) have made important discoveries of transiting low- 
mass planets. Microlensing search es have detected two 
super-Earth s orbiting distant stars (jBeaulieu et al.l[2006l : 
lBennett"et al. 2008) and the statistics of microlensing de- 
tections suggest than cold Neptunes a.re a factor of thre e 
more common that cold Jupiters (jSumi etall [20091) . 
From space, CoRoT has foun d a system with tw o super- 
Earths (one of which transits: iLeger et al.ll2009l ) and Ke- 
pler is poised to detect true Earth analogues in 1 AU or- 
bits using transit photometry with a precision of 20 ppm 
in 6.5 hr (Borucki ct al. 2009). In the next decade, SIM 
Lite (Unw in et al.l 12008') will astrometrically character- 
ize essentially all planets down to Earth mass orbiting 
~100 nearby stars, as well as the more massive planets 
orbiting ■^lOOO stars. 

The Eta-Earth Survey plays a unique role in the study 
of low-mass exoplanets. The population of 230 GKM 
stars in the survey is nearly free of statistical bias since 
the stars were not chosen based on their likelihood of 
harboring a planet, but rather on proximity, brightness, 
and chromospheric activity. Each star is observed at 
least 20 times, insuring a minimum detectability thresh- 
old. Thus, the distributions of planet detections and 
non-detections from the Eta-Earth Survey will yield a 
wealth of information about the efficiency and mecha- 
nisms of planet formation as well as the range of subse- 
quent dynamical histories. The 20 survey observations 
per Eta-Earth Survey target are nearly complete and we 
are aggressively re-observing several promising candidate 
low-mass planets. 

Current theor i es of p lanet formation ()Ida fc LinI 120041 : 
iMordasini et al.l l2009f ) are consistent with the mea- 
sured distributions of massive planets (Saturn mass and 
above) , but their predictions for the abundance and prop- 
erties of low-mass planets are only now being observa- 
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TABLE 1 
Stellar Properties of HD 156668 



Parameter 



Value 



Spectral type 

Mv 

B-V 

V 

J 

H 

K 

Distance (pc) 
[Fe/H] 
Teff (K) 
vsini (kms ^ 
logs 

(Lq) 

(Mq) 
R* (Rq) 

Age (Gyr) 
logR^jK 

Prot (days) 
CTphot (mag) 



K3V 
6.480 
1.015 
8.424 
6.593 
6.117 
6.004 
24.5 ±0.5 
+0.05 ±0.06 
4850 ± 88 
0.50±1.0 
4.598 ±0.12 
0.230 ±0.018 
0.772 ±0.020 
0.720 ±0.013 
8.6 ±4.8 
-4.98 
0.23 
51.5 
< 0.002 



tionally tested. In particular, they predict a dearth of 
planets below roughly Saturn mass in orbits inside the 
ice line. Such planets are predicted to be rare because 
once a planet core grows by planetesimal accretion to a 
of threshold mass of ^lOM^, it undergoes runaway gas 
accretion and becomes a gas giant. While these theories 
consistently predict such a "planet desert" , they differ in 
the contours of the desert in M and a as well as whether 
the planets below the critical core mass survive Type I 
inward migration in significant numbers. 

In this context we announce the discovery of 
HD 156668 b, a super-Earth planet with minimum mass 
Mp s'mi = 4.15 A/® and an orbital period of P = 4.6455 d. 
This is the second super-Earth (Msini<lOM0) to 
emerge from Keck observations explicitly for the 
NASA-UC Eta-Earth Survey, the first being HD 7924 b 
([Howard et al.l[2"009b[ ). The remainder of this paper is 
structured as follows. We describe the host star prop- 
erties in §[21 The spectroscopic observations and their 
Doppler reduction are described in §[3] In §[4l we describe 
the detection of the P = 4.6455 d orbit of HD 156668 b, 
and the high-pass filtering of the RV data that was nec- 
essary to obtain good estimates of the orbital parame- 
ters. In §[5] we carefully consider the null hypothesis — 
the non-existence of HD 156668 b — using S'hk measure- 
ments, photometric observations, and FAP analyses. We 
summarize and discuss the implications of this discovery 
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2. STELLAR PROPERTIES 



HP 156668 (HIP 84607) is a K3 dwarf ([Grav et al.l 
|2003() whose properti es are summarized in Table [T] It 
is nearby {d ~ 24.5 pc; 'van Leeuwen"2007l) and relatively 
bright (V^ = 8.424; H0g et al. 2000). With Afy= 6.480 
and B — V = 1.015, the star lies n ear the Hipparcos av- 
erage main sequence as defined bv IWright (2005). 

Using the SME (S pectroscopy Made Easy) LTE spec- 
tral synthesis code ([Valenti fc Fischeil 120051) , we ana- 
lyzed two high-resolution, iodine-free Keck-HIRES spec- 
tra of HD 156668 and found the effective temperature, 
surface gravity, projected rotational velocity, and iron 



abundance ratio listed in Table [T] The errors on these 
four quantities have been conservatively doubled from 
the formal SME parameter uncertainties because the 
stellar atmosphere models are less certai n at low Teg 
(the SME catalog fValenti fc Fischer' 120051) models stars 
down to Tcs — 4800 K). We refined the above param- 
eters and derived stellar mass, radius, and luminosity 
from SME and interpol ated Yonsei-Yale (Y^) isochrones 
([Demargue et al.l |2004[ ) using an iterative process that 
self-consistently ties together the SME and Y^ values of 
logg ([Valcnti ct al. 2009). HD 156668 appears to be a 
typical middle-aged K dwarf. Its slightly super-solar iron 
abundance of [Fe/H] = -1-0.05 ± 0.06 is consistent with its 
location near the average main sequence. 

Measurements of the cores of the Ca II H & K lines 
show that HD 156668 has modest chromospheric activity 
(Fig. [T]). We measured the chromospheric activity in- 
dices S'hk = 0.23 and log R'^u = —4.9 8 using the method 
descri bed IWright et al.l ([2004D and llsaacson fc Fischeil 
([20091 ). The full set of S'hk measurements for all obser- 
vations of HD 156668 does not show a periodicity near 
the planet's orbital of 4.6455 d (see ii l5.3l for additional 
discussion). 

We estimate a rotati on period Pj-ot ~ 4 8 d using P^j^ 
and B — V calibration ([Noves et al.lll984 ). which is con- 



sistent with the 51.5 d peri od measured by APT photom- 
etry (see ^ 15. 2p . Following IWrighd ([20051 ). and based on 
the values of S'hk, My, and B — V , vfe estimate an RV 
jitter of 1.5ms~^. This empirical estimate is based on 
an ensemble of stars with similar characteristics and ac- 
counts for RV variability due to rotational modulation 
of stellar surface features, stellar pulsation, undetected 
planets, an d uncorrected sy s tematic error s in the velocity 
reduction (iSaar et al.lll998t lWrightil2005D . As explained 
in §[31 jitter is added in quadrature to the RV measure- 
ment uncertainties for Keplerian fitting. 

HD 156668 has several important characteristics that 
make it a nearly ideal RV target star. Like other old, 
chromospherically quiet stars with spectral types from 
late G to early K, HD 156668 appears to be near the min- 
imum of astrophysical jitter arising from acoustic oscil- 
lations, granulation, and photospheric activity (see S I5.3I 
and Figure [7]). The star is relatively bright, yielding high 
signal-to-noise spectra in ~4-5 minutes per observation. 

3. HIRES OBSERVATIONS AND DOPPLER REDUCTION 

We obs erved HD 156668 with the HIRES echelle spec- 
trometer ([Vogt et al.lll994D on the 10-m Keck I telescope. 
The 107 observations span five years (2005-2009) with 
high-cadence observations — clusters of observations on 
6-12 consecutive nights — beginning in 2007. All obser- 
vations were made with an iodine cell mounted directly 
in front of the spectrometer entrance slit. The dense 
set of molecular absorption lines imprinted on the stel- 
lar spectra provide a robust wavelength fiducial against 
which Doppler shifts are measured, as well as strong con- 
straints on the shape of the spectromet er instrumental 
profile at t he time of eac h observation (fMarcv fc Butleil 
[i9_92; Val enti et al. l [T995h . 

We measured the Doppler shift from each star-times- 
iodine spectrum using a mod eling procedure mo dified 
from the method described bv iButler et al.l ([1996D . The 
most significant modification is the way we model the in- 
trinsic stellar spectrum, which serves as a reference point 



A Super-Earth Orbiting HD 156668 



3 




3960 3965 3970 
Wavelength (A) 



3975 



Fig. 1.— Ca II H line for HD 156668 from a Kcck/HIRES spec- 
trum. The line core emission near 3968 A indicates modest chro- 
mospheric activity. 



for the relative Doppler shift measurements for each ob- 
servat ion. Butler et aL use a version of the iJanssonI 
()1995[ ) deconvolution algorithm to remove the spectrom- 
eter's instrumental profile from an iodine-free template 
spectrum. We instead use a new deconvolution algo- 
rithm developed by one of us (J.A.J.) that employs a 
more effective regularization scheme, which results in sig- 
nificantly less noise amplification and improved Doppler 
precision. 

Figure [5] shows RV time series for four stable stars 
with characteristics similar to HD 156668, demonstrating 
our measurement precision of ~1.5-2.0ms~^ (including 
astrophysical, instrumental/systematic, and photon- 
limited errors) for chromospherically quiet late G/early 
K stars over the past 5 years. All of the measurements 
reported here were made after the HIRES CCD upgrade 
in 2004 August and do not suffer from the higher noise 
and systematic errors that limited the precision of pre- 
upgrade measurements to ~2-3ms~^ for most stars. 

The velocities derived from the 107 observations have 
a median single measurement uncertainty of 1.01 ms~^. 
This uncertainty is the weighted standard deviation of 
the mean of the velocity measured from each of the 
^700 2 A chunks in each echelle spectrum (jButler et al.l 
[rgge"). in a few cases, we made consecutive observations 
of HD 156668 to reduce the Poisson noise from photon 
statistics. For the Keplerian orbital analysis below (§|1]), 
the velocities were binned in 2 hr intervals, yielding 86 
measurements with an rms of 2.71 ms^^ about the mean 
and a median measurement uncertainty of 0.99 ms~^. 

4. ORBITAL ANALYSIS 

The measured radial velocities of HD 156668 are listed 
in Table |3] and plotted as a time series in Figure [3^. Be- 
fore any fitting or filtering, these velocities already have 
a low rms, (T = 2.71ms~^. However, long-term trends, 
coherent over at least several months, are apparent by 
visual inspection of Figure [3^ (see, e.g., the downward 
trend in 2007). These tren ds manifest themselves in a 
Lomb-Scargle periodogram (lLomblll976t IScarglel[T98l) of 
the data as long-period power concentrated near ~2.3yr 
(Figure [3]3). 

The RVs also show a significant periodicity near 4.647 d 
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Fig. 2. — Radial velocity time series for four stable stars in our 
Keck Doppler survey that are similar to HD 156668. These stars 
demonstrate long-term velocity stability and precision of HIRES. 
The binned velocities with measurement uncertainties (but not jit- 
ter) are plotted. Panels are labeled with star name, spectral type, 
and velocity rms. 

with power than exceeds the 0.1% FAP threshold^°. 
A second short-period peak is visible at 1.270 d, but 
this peak is a stroboscopic alias of the 4.647 d sig- 
nal with the sidereal day length (see ? I5.4I) . We at- 
tempted to fit the radial velocities with a single-planet 
Keplerian orbital solution seeded with P = 4.647 d us- 
ing a partially-linearized, least-squares fitting procedure 
(jWright & Howard 2009). Each velocity measurement 
was assigned a weight constructed from the quadra- 
ture sum of the measurement uncertainty (listed in 
Table [3]) and a stellar jitter term (1.5 ms^^). The 
fitting routine converged on a robust solution with 
P = 4.6455 d, e = 0.39, and ii: = 2.28ms-^ This model 
gives a = 2.33 ms~^ and Xi^ = 1-70, a significant improve- 
ment over a linear model to the data. 

We also tried fitting the long-period signal with a 
single-planet Keplerian model. We seeded the fitting rou- 
tine with a range of periods and found a best-fit solution 
with P = 791 d (2.17yr) and a poorly constrained eccen- 
tricity when that parameter is allowed to fioat. The im- 

The false alarm probability (FAP) thresholds plotted as 
dashed horizontal lines in the periodograms in Figure [3] refer to 
the probability that periodograms of random rearrangements of the 
data would exceed the specified power level. Since periodograms 
only measure the power of sine wave fits to the data (i.e. circular 
orbits), these FAPS are less conservative than the ones described 
in ij lS.ll that allow for eccentric orbits and use the AxJ statistic. 
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Fig. 3. — RV measurements and 5hk values for HD 156668 based on Keck-HIRES spectra. Each pair of panels in a given row displays 
the same data as a time series (panels a-d, on left) and as a periodogram (panels e-g, on right). Panels a and e show the unfiltered RVs 
from Table [3I with significant Fourier power at ~2.i5yr, 4.647 d, and 1.270 d (an alias of 4.647 d with the sidereal day length). Panels b and 
/ show the RVs after applying a high-pass filter by subtracting the dashed line model in panel a; note that the ~2.3yr signal was almost 
completely removed, strengthening the power at 4.647 d and its alias. Panels c and g show the residuals of the high-pass filtered data to 
a one-planet fit with P = 4.6455 d in a circular orbit. Panels d and h show values of the chromospheric index Shk derived from each RV 
measurement. 



provements in a and Xiy for this model are comparable 
to those from the P = 4.6455 d model described above. 
However, the poorly constrained and high eccentricity of 
the fit and the uneven phase coverage render the inter- 
pretation of this signal suspect. We cannot be sure if it 
is due to a planet in Keplerian motion, an astrophysical 
signal such as chromospheric activity masquerading as a 
RV change of the star, or some other effect. 

Whatever the source of the long-period signal, its ef- 
fect is to obscure the short-period signal. We consid- 
ered several ways to high-pass filter the RVs to isolate 
the short-period signal. A key requirement of such a 
filtering process is that it leave the short-period sig- 
nal untouched. We concluded that fitting the data 
to a two-planet model plus a linear trend would ro- 
bustly separate the long and short-period signals, al- 
lowing for accurate parameter estimation for the short- 
period signal and a fair assessment of its statistical 
significance. The fitting routine converged on a ro- 
bust solution with Keplerian parameters for the inner 
planet listed in Table [21 We adopt these parameters 
for HD 156668 b. The "outer planet" in this model 



TABLE 2 
Orbital Solution for, HD 156668 b 



Parameter 



Value 



P (days) 

Tc (JD - 2,440,000) 

e 

K (ms-i) 
Msini (M®) 
a (AU) 
^oba (binned) 

Median binned uncertainty (ms" 
Assumed jitter (ms~^) 
(T (ms"-"-) 



4.6455 ± O.OOll 
14718.57 ± 0.11 
=0.0 
1.89 ± 0.26 
4.15 ±0.58 
0.0500 ±0.0007 
86 
0.99 
1.50 
1.74 
0.97 



has P = 2.31yr, e = 0.48, if = 4.05 m s" 



U! - 



178 deg. 



Tp^ 2,454,414.22, and dw/di = 0.94m s'^ yr"^ as shown 
by the dashed line in Figure [3l We emphasize that the 
outer planet in this fit is only a convenient model for the 
data, effectively serving as a high-pass filter to isolate the 
signal of the inner planet. Determining whether the long- 
period signal represents a planet will require additional 
RV measurements and diagnostic data. 
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Figure Eb/f shows the RV measurements after sub- 
tracting the long-period signal. Trends in the time se- 
ries are no longer apparent and the 4.647 d signal and 
its ahas at 1.270 d are significantly strengthened in the 
periodogram. Figure [3J;/g shows the RVs after subtract- 
ing both long and short period signals. The value of a is 
1.74ms~^ and the periodogram appears nearly feature- 
less, with no periods remaining having significant power. 

The phased orbital solution for HD 156668 b is 
shown in Figure SI The Doppler semi-amplitude of 
X = 1.89 m s""'^ is extremely low and is nearly equal 
to the 1.80 ms~^ typical error for single measurements 
(including jitter). The resulting minimum mass of 
M sin i = 4.15 is also extremely small, and is the sec- 
ond lowest reported to date using the RV technique. We 
adopted a circular orbit because the uncertainty on the 
eccentricity (e = 0.20 ± 0.17) is significant when that pa- 
rameter is allowed to float in the two-planet fit. Further, 
the best-fit eccentric orbit model did not show a statis- 
tically significant improvement in compared to the 
circular orbit model. 

The Keplerian parameter u ncertainties were derived 
using a Monte Carlo method (|Marcv et al. 2005). The 
uncertainty estimates on M sin i and a account for the 
uncertainty from M^, . With five years of observations of 
this short-period planet, the error on P is quite small 
(one part in 4000). 

5. THE NULL HYPOTHESIS 

In this section we explicitly consider the null hypothe- 
sis. That is, we consider whether the P = 4.6455 d signal 
we identified as due to a low-mass planet in Keplerian 
motion could be due to another source. Such an anal- 
ysis is motivated by the extremely low amplitude signal 
(ivT = 1.89 m s~^) and because of the use of high-pass fil- 
tering. 

5.1. False Alarm Probabilities 

We considered the possibility that the P = 4.6455 d sig- 
nal arose from the chance arrangement of random, statis- 
tically independent errors i n the RVs by computing FAPs 
for several fits to th e data f Marcv et al.ll20d5t iCummiii^ 
|2004; Ho ward et al. 2009a^b). These FAPs compare the 
measured data to 1000 scrambled data sets drawn ran- 
domly with replacement from unscrambled data. For 
each data set we compare a best-fit Keplerian model to 
the null hypothesis (a linear fit to the data) by comput- 
ing Axl = xL,u - Xkcp,u' wliere xL,^ and xhp,u are the 
values of xt for linear and Keplerian fits to the data, re- 
spectively. The statistic measures the improvement 
in the fit of a Keplerian model compared to a linear model 
of the same data. The FAF is the fraction of scrambled 
data sets that have a larger value of Axi than for the 
unscrambled data set. That is, the FAP measures the 
fraction of scrambled data sets where the improvement 
in Axi from a best-fit Keplerian model over a linear 
model is greater than the improvement of a Keplerian 
model over a linear model for the actual measured veloc- 
ities. We use Axl as the goodness-of-fit statistic, instead 
of other measures such as Xi' ^or a Keplerian fit, to ac- 
count for the fact that the scrambled data sets, drawn 
from the original velocities with replacement, have differ- 
ent variances, which sometimes artificially improve the 



fit quality (i.e. some scrambled data sets contain fewer 
outlier velocities and have lower rms). Note that this 
FAP does not measure the probability of non-planetary 
sources of true velocity variation masquerading as a plan- 
etary signature. 

We computed an FAP for the "two-planet model" 
(with the long-period planet acting as a high-pass fil- 
ter) by comparing Axl for the actual data with Axl for 
two-planet fits to scrambled sets of the high-pass filtered 
data with the unscrambled long-period signal added back 
in. This FAP has a value of 0.3% (3/1000) and tests 
whether the statistical significance of the short-period 
signal is somehow enhanced by simultaneously fitting for 
a long-period signal. 

We also computed the FAP for a 1-planet fit to the 
high-pass filtered data. We compared Axl for the 
unscrambled, filtered data set and scrambled versions 
thereof. This FAP tests whether the short-period signal 
is a statistical fluctuation, while assuming that the long- 
period signal is real (but still of unknown origin). We 
found an FAP of 0.4% (4/1000) for this scenario (with 
eccentricity unrestricted in the fits to scrambled data 
sets). We note that all four of the false alarm solutions 
had best-fit parameters that appeared unphysical: high 
eccentricity (e>0.5) and short period (P<10d). (We 
regard them as unphysical because short-period planets 
are almost universally in near circular orbits.) Such spu- 
rious, high eccentricity solutions often appear as the best- 
fit solution to low rms RV data lacking a coherent signal. 
When we restricted the fits of scrambled data sets to cir- 
cular orbits (as we did with the actual measurements), 
the FAP dropped to <0.1% (0/1000). 

We conclude that the P = 4.6455 d signal is statistically 
significant. Considering all of the FAP tests, we estimate 
that the probability that this signal arose just due to 
random errors is less than 1% and probably on the order 
of 0.1%. 

5.2. Photometric Observations 

In addition to the Keck RVs, we obtained contempora- 
neous Stromgren b and y photometric measurements with 
the TIO 0.80 m automatic photometric telescope (APT) 
at Fairborn Observatory in Arizona. This APT is identi- 
cal to the T8 APT that was used to acquire photometric 
observations of HD 79 24 in the first paper of this series 
([Howard et al.l[2"009b[ ). (That paper mistakenly identi- 
fied the APT used as the T12 0.8 m APT.) Our observing 
procedures and data reduction p rocedures in this pape r 
are identical to those described in lHoward et al.l (|200 9bD. 

Our three comparison stars A, B, and C for HD 156668 
(star D) were HD 158974 {V = 5.63, B -V = 0.96, 
G8 HI), HD 155092 {V = 7.07, B -V ^ 0.42, F2), and 
HD 156536 {V = 7.51, B - V = 0.42, F3 IV), respec- 
tively. The TIO APT acquired 477 good observations of 
this quartet of stars covering the 2007, 2008, and 2009 ob- 
serving seasons. Comparison star A, HD 158974, exhibits 
low-amplitude brightness variability of 5 mmag with a 
period of 152 days. The C — B differential magnitudes 
have a standard deviation of 1.7 mmag, indicating that 
both stars are constant to the level of our measurement 
precision. Therefore, we created differential magnitudes 
of HD 156668 by averaging the D — B and D — C differ- 
ential magnitudes into a single D — {B + C)/2 differential 
magnitude to be used for our analyses. We also combined 
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Fig. 4. — Single-planet model for the radial velocities of HD 156668, as measured by Kcck-HIRES. The dashed line shows the best-fit 
circular orbital solution representing a 4.15 M^g (minimum mass) planet in a 4.6455 d orbit. Filled circles represent phased, binned, and 
high-pass filtered velocities, while the open circles represent the same velocities wrapped one orbital phase. The error bars show the 
quadrature sum of measurement uncertainties and 1.5 ms^^ jitter. 



our b and y observations into a single {b + y)/2 passband 
to improve the precision further. 

These D — [B + C)/2 magnitudes in the {b + y)/2 pass- 
band are shown in the top panel of Figure [5] The gaps 
following the longer runs of data are due to the neces- 
sity of shutting down APT operations during the summer 
rainy season in Arizona. The gaps following the shorter 
data groups are the normal seasonal gaps for HD 156668. 
A summary of the photometric observations is given in 
TableO Column 5 indicates that the seasonal mean mag- 
nitudes vary over a range of '^Immag, which is typical 
for solar- type stars with modes t chromospheric activity 
(see, e.g.. iLockwood et al.ll2007D . 

The vertical bars in the top panel of Figure [5] encom- 
pass the portion of the D — {B + C)/2 light curve that 
most clearly exhibits coherent variability that might be 
due to rotational m odulation of spots on the star's pho- 
tosphere (see, e.g.. iHenrv et al.l llOQSl. That section of 
the light curve is replotted in the second panel. Straight 
line segments approximate the 2.5 cycles of the purported 
brightness variability. 

Plotted in the third panel of Figure [5] is the frequency 
spectrum of the observations from panel two, showing a 
clear periodicity of 51.5 d. In the bottom panel, the same 
observations are phased to the 51.5 d period and a time 



of minimum computed with a least-squares sine fit to 
the observations. The sine fit also gives a peak-to-peak 
amplitude of 2.3 ± O.Smmag. 

We take the 51.5 d period to be the rotation period 
of HD 156668. This is consistent with v sin z = 0.5 km 
and is very close to the rotation period of 48 d predicted 
from its level of chromospheric activity. 

The APT photometry is also useful for confirming that 
observed RV variations are due to a planetary compan - 
ion and not stellar surfac e activity. iQueloz et al.l (|2001l ) 
and 'Pauls on et al.l (|2005f ) have demonstrated how rota- 
tional modulation in the visibility of star spots on ac- 
tive stars can result in periodic RV variations and their 
misinterpretation. All 477 photometric observations of 
HD 156668 are plotted in the top panel of Figure [SI 
phased with the 4.6455 d RV period and a time of mid- 
transit computed from the orbital elements. A least- 
squares sine fit on that period gives a semi-amplitude of 
only 0.10 ± O.llmmag. This absence of detectable ro- 
tational modulation of surface activity to high precision 
on the RV period provides strong evidence that the RV 
variations arise from a super-Earth companion. 

The bottom panel of Figure [5] replots the photometric 
observations that lie near the predicted time of transit. 
The solid curve shows the predicted depth (S.lmmag) 
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TABLE 1 

Summary of Photometric Observations of HD 156668 



Observing 




Date Range 


Sigma 


Seasonal Mean 


Season 


Nobs 


(HJD - 2,400,000) 


(mag) 


(mag) 


(1) 


(2) 


(3) 


(4) 


(5) 


2007 


181 


54117-54407 


0.00195 


1.27046 ±0.00014 


2008 


225 


54484-54754 


0.00188 


1.27121 ±0.00013 


2009 


71 


54865-55004 


0.00131 


1.27013 ±0.00016 



1.27 



+ 

<" 1.275 




54500 54550 54600 

JD - 2,400,000 



3 
t3 



0.3 




t = 0.01942 


0.2 




P = 51.5 : 


0.1 














0.1 0.2 0.3 

Trial Frequency (cycles/day) 




Photometric Phase 

Fig. 5. — Top: The 477 D-{B + C)/2 photometric observations 
of HD 156668 in the {b + y)/2 passband, acquired with the TIO 
0.8 m APT during the 2007, 2008, and 2009 observing seasons. 
Second panel: The portion of the 2008 observing season set off 
by vertical bars in the top panel shows the most coherent bright- 
ness variability due to cool star spots carried across the disk of 
the star by its rotation. Third panel: Frequency spectrum of the 
observations in the second panel gives a best period of 51.5 days. 
Bottom Panel: Plot of the data from panel two phased with the 
51.5-day period reveals coherent varialaility with a peak-to-peak 
brightness amplitude 0.0023 mag. 



and duration (±0.011 phase units) of a central transit, 
computed from the orbital elements and stellar proper- 
ties with an assumed planetary composition of pure hy- 
drogen. That such shallow tran sits can be detecte d with 
the APTs has been shown by iSato et all ()2005[) . who 
used the Til APT to discover the 3mmag transits of 
HD 149026b. The dotted line across the transit window 
corresponds to the expected depth (0.6 mag) of a hypo- 
thetical planet composed of entirely water. Our observa- 
tions show that transits with a depth of 3.1 mag probably 
due not occur, essentially ruling out transits of a hydro- 
gen planet. Additional precise photometry is required to 
rule out other planetary compositions. 



C\j 1.265 



1.275 
1.28 



0.8 



0.0 



0.2 



0.4 



0.6 



0.8 




0.94 



0.96 



0.98 



0.00 



0.02 



0.04 



0.06 



Planetary Orbital Phase 



Fig. 6. — Top panel: The 477 APT photometric observations of 
HD 156668 plotted modulo the 4.6455 d period of the RV variations. 
A least-squares sine fit at that period yields a semi-amplitude of 
only 0.10 ± 0.11 mmag, providing strong evidence that the velocity 
variations arise from a planetary companion. Bottom panel: The 
photometric observations of HD 156668 near the predicted time of 
transit replotted with an expanded scale. The solid curve shows 
the predicted depth (3.1 mmag) and duration (±0.011 phase units) 
of a central transit, computed from the orbital elements, the stel- 
lar properties, and a planetary composition of pure hydrogen. The 
uncertainty in the predicted transit time is shown by the error 
bar under the transit window. The dotted line across the transit 
window shows the expected depth (0.6 mmag) for a planet com- 
posed of water. Transits of a pure hydrogen planet are essentially 
ruled out. Additional observations are required to rule out other 
planetary compositions. 



5.3. Chromospheric Activity 

RV planet searches measure the shifts of the centroids 
of thousands of stellar lines. The shapes of individual 
lines are determined in large part by Doppler broaden- 
ing from the disk-averaged velocity field of the star's 
surface. Spots are magnetically-controlled regions of 
the stellar photosphere characterized by temperatures 
~1000K lower than unaffected regions. As they rotate 
across the stellar surface, spots contribute less flux to 
particular parts of each absorption line profile — less flux 
on the blue side of lines for spots on the approaching 
limb and less flux on the red side for spots on the re- 
ceding limb — thereby distorting the average line profile, 
shifting the centroid, and causing an apparent Doppler 
shift of the star. 

False positive signals of this type te nd to occur 
around chromospherically active stars fOueloz et al.l 
i2001t iPaulson et al..,2005i ). The spurious RV signals are 
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Fig. 7. — Plot of logij^p, as a function of B — V color for all G 
and K stars in the Eta-Earth survey with B-V <1.2. HD 156668 
(large red triangle highlighted by dashed red lines) is among the 
least chromospherically active Eta-Earth stars with similar B — V. 

coherent over typical spot lifetimes (weeks to months) 
and have periods similar to the stellar rotation pe- 
riod. In S I5.2I we showed that the 4.6455 d signal is not 
seen in APT photometry, ruling out rotationally mod- 
ulated spots as the source of the RV periodicity. Here, 
we use the chromospheric indices S'hk and log -R^k 
strengthen that conclusion further. These indices mea- 
sure the level of stellar chromospheric activity, which in 
turn is strongly correlated with the magnetic activity of 
the stellar photosphere. 

Some stars also show long-term chromospheric activ- 
ity cycles as the average number of spots rises and falls 
with the solar cycle, typically with a timescale of ~10yr. 
These cycles are sometimes detected as apparent RV 
shifts and incorrectly interpreted as long-period planets. 

Figure [7] shows the average value of logi?^j^ for 
HD 156668 and the other stars in the Eta-Earth Survey 
with B — V < 1.2. HD 156668 is among the least chro- 
mospherically active stars of similar spectral type in the 
Eta-Earth Survey. This is consistent with our limits on 
photometric variability from APT measurements fii l5.2l) 
and the small RV residuals to the Keplerian fit. 

The values of S'hk derived from each Keck-HIRES 
spectrum are listed in Table [3] and plotted as a time 
series in Figure [3ji. A Lomb-Scargle periodogram of 
the data is plotted in Figure [3]i. Importantly, the pe- 
riodogram shows negligible power for P < 20 d, as ex- 
pected for a middle-aged star with a predicted rotation 
period of ^48 d. This lack of power at periods near 4.6 d 
strengthens the case for HD 156668 b. 

On average, S'hk rises with time (Figure [3ji) , possibly 
due to the observation of a partial solar magnetic cycle. 
This several year long trend is not observed in the RVs 
(Figure [3^) . One feature of the RV time series is also 
apparent in the Shk, the declining trend in late 2007. 
This apparent correlation raises the possibility that the 
long-period signal that was filtered out of the RVs in §2] 
is due to stellar activity. Despite this apparent correla- 
tion, the Pearson linear correlation coefficient, r — -t-0.11, 
between the unfiltered RVs and the Shk values indicates 
an insignificant correlation. When we subtract a second- 
order polynomial fit to the Shk values, the correlation 
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Fig. 8. — Running periodogram power near P = 4.6455 d as 
a function of time. Each dot represents the maximum peri- 
odogram power for the set of all high-pass filtered RVs up to 
that point in time. Only periods near the orbital period of 
HD 156668 b (P = 4.64-4.65 d) were used to compute each maxi- 
mum periodogram power value. The monotonic rise confirms that 
the signal was present throughout the observations and points to 
a dynamical origin. 

with the RVs is still statistically insignificant, r = -t-0.19. 
Thus, the suggestion that stellar activity explains the 
long-period signal is not well supported when the entire 
data set is considered. 

Rotational modulation of spots and other surface fea- 
tures can also be d etected in the spectral line bisectors 
(|Torres et al.ll2005l) . However, for this extremely low am- 
plitude signal {K <2ms~^), detecting line profile varia- 
tions at the same or higher precision is not possible. 

As a final check, we show in Figure |8] that the RV 
signal from HD 156668 b is strictly periodic and present 
throughout the observations, as the clock-like signal from 
a planet should be. The plot shows the periodogram 
power of the planet rising monotonically as additional 
measurements are taken. False periodicities, such as 
those due to spots, typically exhibit periodicities that 
are only briefly coherent and that may reappear with 
slightly different periods. 

5.4. Stroboscopic Alias 

We interpret the periodogram spike near 1.270 d (Fig- 
ure [3]) as a stroboscopic alias of the orbital period (which 
appears on the periodogram as 4.647 d due to finite pe- 
riod sampling) with the sidereal day length (0.997 d): 
l/4.647d + l/1.270d = l/0.997d. We verified this by 
computing the periodograms of 2000 mock data sets that 
have the same times of observations and errors as the ac- 
tual RVs. Of the 1000 data sets that we injected a 4.647 d 
signal, 851 of them showed significant periodogram power 
(> 10) near 1.27d. Conversely, 0/1000 of the mock data 
sets with a 1.270 d signal added showed significant power 
near 4.647 d. 

6. DISCUSSION 

We present the detection of HD 156668b, a super- 
Earth planet with minimum mass Mp sin i = 4.15 
in a P = 4.6455 d orbit around a K3 dwarf. We draw 
on several lines of evidence to support the existence of 
HD 156668b. We showed in §[5T] that the short-period 
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signal is statistically significant. This signal is apparent 
in a fit to the unfiltered RVs, and stands out strongly 
when isolated by high-pass filtering. The host star is 
middle-aged and quiet, providing a nearly ideal RV tar- 
get. The planet's short-period signature is not seen in 
photometric observations or in chromospheric indices. 
Thus, the evidence strongly points to a planetary inter- 
pretation. 

To estimate the orbital parameters of HD 156668 b, 
we found it necessary to apply a high-pass filter to the 
RVs by subtracting the model of a long-period Kcple- 
rian. This model changes slowly over time and clarifies, 
but does not artificially enhance, the i-' = 4.6455 d sig- 
nal of HD 156668 b. Filtering of this type is common 
in other areas of time domain astronomy (e.g. transit 
photometry, astroseismology) and is uncommon, but not 
unprecedented, in RV planet detection. Similar filtering 
techniques were em ployed to disen t angle the planetary 
signals of G1 176b (iF^ eille et all [2009h and G1674b 
(|Bonfils et al.l 120071 ) . although in those cases the non- 
planetary signal wa s clearly due to sp ots modulated 
by stellar rotation. iQueloz et al.l (|2009l ) used a "pre- 
whitening" technique to extract the signatures of two 
super-Earths from the complicated RV time series of the 
significantly more active star, Corot-7. Further, many 
exoplanets are announced with a model that includes a 
linear or second order RV trend, presumably due to a 
long period orbital companion. The source of the long- 
period signal remains unknown and motivates additional 
measurements. If it is due to a planet, the body has 
minimum mass ^^45 Mq and orbits with Pc — 2.31 yr and 
ac = 1.6 AU, a cold super-Neptune near the ice line. 

We see no evidence for transits of HD 156668 b down 
to a level of ~3mmag. However, given the large a priori 
transit probability of 7%, it is instructive to speculate 
about the transit signatures of va rious possible plane t 
compositions. Using the models in iSeager et al.l ()2007[ ). 
a 4M0 planet composed of pure Fe, MgSiOa, H2O, or 
H would yield planets of radius i?pi = 1.2, 1.5, 2.0, and 
4.5 i?®, producing transits of depth 0.22, 0.35, 0.61, and 
3.1 mmag, respectively. These homogeneous planet mod- 
els are oversimplified, but set the scale for admixtures 
of those ingredients: transits of planets made of solids 
and water would have depths of ~0. 2-0. 6 mmag, while 
transits of a planet with a significant atmosphere could 
be much deeper. For comparison, the transiting super- 
Earth Corot-7 b has a transit depth of 0.36 mm ag im- 
plying a radius i?pi = 1.68 i?© (iLfeer et al I I2009D. Us- 
ing the RV-determined mass Mpi = 4.8 M^p ()Queloz et ahl 
I2009D . the bulk density is terrestrial, ppi = 5.6gcm~^. In 
contrast, GJ 1214b has Afpi = 6.6Me and i?pi = 2.7i?©, 
implying ppi = 1.9gcm~'^, intermediate between Earth 
and the ice giants of the Solar System. Transits of 
GJ 1214 b are unusually deep (15 mmag) for a planet of 
this size because it orbits an M dwarf with i?^, = 0.21 Rq 
(jCharbonneau et al.ll2009l). 

Sev eral authors fe.g. Ilda fc Lin|[200^ iMordasini et al.l 
|2009| ) have argued that super-Earths will have insignif- 
icant hydrogen atmospheres (by mass) because during 
formation their masses failed to reach a critical mass 
(typically ^lOAf©) when the growth from solid plan- 
etesimals is augmented by runaway accretion of gas from 
the protoplanetary disk. Smaller atmospheres (up to 



several per cent by mass) could be produced by de- 
gassing during impact accretion and geological activity 
(jElkins-Tanton fc Seageil [2008t iKitreTaLl l2009l) . How- 
ever, whatever atmosphere is acquired fr om these pro- 
cesses ma y be lost t o atmospheric escape (jBaraffe et al.l 
2006.; Valencia et al.. 2009). Nevertheless, the brief his- 
tory of exoplanets is replete with observational sur- 
prises (hot Jupiters, eccentric orbits, etc.) so we con- 
sider the obs e rvatio nal consequences of an atmosphere. 
lAdams et al.l (|2008| ) find that adding a H/He gas enve- 
lope equivalent to 0.2-20% of the mass of a solid 5 M© ex- 
oplanet increases the radius 8-110% above the gas-free 
value. Atmospheres dominated by heavier molecules 
such as H2O and N2 (as on Earth) would swell the planet 
less for the same atmospheric mass because of the higher 
mean molecular weight and reduced scale height. Thus, 
we conclude that the APT photometric observations rule 
out transits for HD 156668 b if the radius is dominated 
by a H/He atmosphere (tens of per cent by mass), but 
do not preclude transits if the atmosphere is less massive 
or composed of heavier elements. 

HD 156668 b adds statistical weight to the emerging 
trends of the properties of super-Earths and their hosts. 
Like most other such planets, it orbits a K or M dwarf. 
In contrast to other super-Earth hosts, HD 156 668 has 
a slightly super-solar metallicity (|Howard et al.l [200 9b). 
The rate of multiplicity in systems with super-Earths 
and Neptune-mass planets appears to be much higher 
than f or higher mass planet hos ts, with H D 40307 
(iMavoret al,, ,2009a). G J581 ((Mavor et al.l[20b9b ). and 
HD 69830 ()Lovis et al.ll2006[ ) being the standard exam- 
ples of multiplicity in low-mass systems. The long-period 
signal seen for HD 156668 is suggestive of a second low- 
mass planet in the system. 

HD 156668 b pushes the frontier of RV planet dis- 
covery to lower masses and smaller Doppler ampli- 
tudes. It is the second lowest minimum mass exo- 
planet discov ered to date by th e RV technique, af- 
ter GJ581e ()Mavor et al.l l20091j ). With a Doppler 
semi-amplitude of 1.89ms~^, HD 156668b is also only 
the second exoplanet discovered to date with K < 
2.0ms"^ GJ 581 e is the other with if = 1.85m s"i and 
Msini = 1.9M(5 (|Mavor et all l2009bD . This progress is 
remarkable: 51Pegb was discovered (jMavor fc Queloj 
I1995D with Af sin 1 = 0.47 Mjup and K = b7-ms-^ while 
the signal from HD 156668 b is smaller by factors of 35 
and 30, respectively. However, we expect that true Earth 
analogs — 1 planets in 1 AU orbits around G stars — 
will remain permanently out of reach of the Doppler tech- 
nique and will instead be discovered by transit photom- 
etry, astrometry, and microlensing. The Earth imparts a 
i<r = 0.09ms~^ signal on the Sun, a factor of 20 smaller 
than HD 156668 b. Our Keck/HIRES observations of 
HD 156668 take ~5min to achieve 1 ms~^ precision; 
scaling to 0.1 ms^^ precision requires 100 times the pho- 
tons and impractically long integration times. (A hand- 
ful of exceptionally bright nearby stars such a s a Cen 
A/B may be an exception ()Guedes et al.ll200"8l) .) More- 
over, stellar jitter, even when averaged over with long in- 
tegrations, likely limits Doppler precision to a few tenths 
ofams""'^ for the most quiet stars. Longer period planets 
suffer from the additional disadvantage that only a few 
orbits transpire during an observational campaign (com- 
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pared with ~100 orbits for the discoveries of HD 155668 b 
and GJ581e), reducing their detectabiUty by clock- Uke 
coherence. 

Nevertheless, the future of the Doppler technique is 
bright and we expect that it will continue to play a 
prominent role in planet searches. Discoveries of super- 
Earths and Ncptuncs from the Eta-Earth Survey will 
shape our understanding of planet formation and mi- 
gration. Dopplcir work in other domains long period 
giant planets, subgiants, multi-planet systems, dynami- 
cally interacting planets, Rossiter-McLaughlin measure- 
ments, etc. — will also continue to inform and enrich our 
knowledge of the rich astrophysics of exoplanets. 
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TABLE 3 

Radial Velocities and Shk values for HD 156668 



Radial Velocity Uncertainty 
JD- 2440000 (ms-i) (ms^i) Shk 



13478.97768 



-1.04 



1.03 



0.202 



A Super-Earth Orbiting HD 156668 



TABLE 3 — Continued 





Rcidicil Velocity 


Uncertciinty 




TD — 9440000 


(ms -*-) 


I'm ii-l'l 




13547.90964 


-4.86 


0.99 


0.214 


13604.83890 


-5.21 


0.85 


0.213 


13807.14411 


3.29 


1.14 


0.222 


13932.91866 


0.87 


1.00 


0.229 


13960.91401 


3.29 


0.87 


0.220 


13961.80956 


0.06 


0.88 


0.219 


13981.77064 


3.87 


0.94 


0.231 


13982.87659 


0.94 


0.84 


0.233 


13983.81911 


4.28 


0.80 


0.230 


13984.90602 


5.74 


0.85 


0.231 


14247.01977 


1.27 


0.71 


0.239 


14248.00711 


1.47 


1.10 


0.238 


14249.92154 


0.79 


1.14 


0.235 


14252.01612 


4.28 


1.08 


0.234 


14255.84743 


-0.27 


0.76 


0.238 


14277.79111 


-3.66 


1.11 


0.230 


14278.80991 


-2.37 


0.98 


0.230 


14285.81562 


0.75 


1.10 


0.230 


14294.89361 


2.08 


1.16 


0.226 


14304.95027 


-0.99 


0.75 


0.232 


14305.95179 


-2.31 


0.69 


0.234 


14306.92739 


-2.31 


0.70 


0.232 


14307.97764 


-0.33 


0.73 


0.231 


14308.94701 


2.07 


1.00 


0.229 


14309.94163 


0.68 


0.71 


0.232 


14310.93496 


0.77 


0.66 


0.230 


14311.92755 


-1.40 


0.67 


0.230 


11312.93048 


1.49 


0.72 


0.228 


14313.92787 


2.74 


0.68 


0.228 


14314.96879 


-2.10 


1.01 


0.227 


14335.87372 


-2.54 


0.63 


0.229 


14343.79865 


-4.20 


0.89 


0.219 


14396.70306 


-3.13 


0.95 


0.221 


14397.70527 


-4.94 


0.93 


0.222 


14398.73491 


-6.92 


1.11 


0.219 


14636.95177 


-3.36 


1.23 


0.238 


14717.74461 


3.34 


0.99 


0.254 


14718.89500 


-0.22 


0.97 


0.252 


14719.78911 


-4.40 


1.02 


0.256 


14720.82821 


2.28 


0.92 


0.255 


14721.81321 


5.53 


1.00 


0.255 


14722.75844 


0.73 


0.97 


0.255 


14723.75113 


-0.66 


0.90 


0.255 


14724.77287 


-2.25 


0.96 


0.254 


14725.74632 


0.03 


1.03 


0.252 


14726.80657 


4.30 


0.97 


0.250 


14727.81096 


1.82 


0.92 


0.250 


14777.68651 


2.56 


1.05 


0.236 


14779.69190 


-4.18 


1.23 


0.235 


14927.05569 


-1.01 


1.31 


0.251 


14930.09564 


2.92 


1.26 


0.244 


14935.05529 


4.30 


1.19 


0.237 


14955.10678 


4.61 


1.16 


0.239 


14955.97159 


0.56 


1.14 


0.241 


14957.05565 


-1.73 


1.05 


0.241 


14983.97707 


-0.72 


1.09 


0.241 


14984.96665 


-2.64 


1.22 


0.242 


14985.92186 


0.69 


1.13 


0.239 


14986.99173 


1.99 


1.14 


0.240 


14987.99574 


-1.07 


1.20 


0.242 


14988.89618 


-2.23 


1.08 


0.244 


15016.97485 


-4.08 


1.08 


0.253 


15019.04148 


0.23 


1.25 


0.252 


15041.97022 


-0.80 


1.16 


0.236 


15042.89409 


2.71 


1.24 


0.235 


15043.90466 


-1.71 


1.13 


0.233 


15044.94394 


-2.75 


1.00 


0.233 


15048.86672 


-2.97 


1.26 


0.239 


15073.76736 


-1.01 


0.95 


0.255 


15074.75962 


1.61 


0.91 


0.253 


15075.73804 


0.05 


0.66 


0.252 


15076.74458 


-2.42 


0.68 


0.247 


15077.74489 


-1.57 


0.96 


0.242 


15078.76561 


-2.64 


1.00 


0.239 


15079.73946 


-2.60 


0,72 


0.237 


15()8l).7i:«)l) 


-0.13 


0.(><) 


0.235 



Howard et al. 



TABLE 3 — Continued 





Rcidicil Velocity 


Uncertciinty 




Tf) _ 244nnnn 


(m s -*-) 


Cm 'i-l'l 




15081.73701 


-3.35 


0.69 


0.234 


15082.72821 


0.83 


1.01 


0.229 


15083.73498 


1.02 


0.69 


0.230 


15084.74428 


1.19 


1.00 


0.228 


15106.76608 


-2.26 


1.16 


0.256 


15109.75425 


-3.25 


1.39 


0.261 


15111.72700 


-1.94 


0.79 


0.263 


15134.69674 


-0.50 


0.82 


0.232 


15135.70259 


-1.59 


0.80 


0.229 



